\(\int \frac {\sec (e+f x)}{(a+b \sec (e+f x)) (c+d \sec (e+f x))} \, dx\) [256]

Optimal result

Integrand size = 31, antiderivative size = 121 \[ \int \frac {\sec (e+f x)}{(a+b \sec (e+f x)) (c+d \sec (e+f x))} \, dx=\frac {2 b \text {arctanh}\left (\frac {\sqrt {a-b} \tan \left (\frac {1}{2} (e+f x)\right )}{\sqrt {a+b}}\right )}{\sqrt {a-b} \sqrt {a+b} (b c-a d) f}-\frac {2 d \text {arctanh}\left (\frac {\sqrt {c-d} \tan \left (\frac {1}{2} (e+f x)\right )}{\sqrt {c+d}}\right )}{\sqrt {c-d} \sqrt {c+d} (b c-a d) f} \]

[Out]

2*b*arctanh((a-b)^(1/2)*tan(1/2*f*x+1/2*e)/(a+b)^(1/2))/(-a*d+b*c)/f/(a-b)^(1/2)/(a+b)^(1/2)-2*d*arctanh((c-d)
^(1/2)*tan(1/2*f*x+1/2*e)/(c+d)^(1/2))/(-a*d+b*c)/f/(c-d)^(1/2)/(c+d)^(1/2)

Rubi [A] (verified)

Time = 0.33 (sec) , antiderivative size = 121, normalized size of antiderivative = 1.00, number of steps used = 6, number of rules used = 4, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.129, Rules used = {4073, 3080, 2738, 214} \[ \int \frac {\sec (e+f x)}{(a+b \sec (e+f x)) (c+d \sec (e+f x))} \, dx=\frac {2 b \text {arctanh}\left (\frac {\sqrt {a-b} \tan \left (\frac {1}{2} (e+f x)\right )}{\sqrt {a+b}}\right )}{f \sqrt {a-b} \sqrt {a+b} (b c-a d)}-\frac {2 d \text {arctanh}\left (\frac {\sqrt {c-d} \tan \left (\frac {1}{2} (e+f x)\right )}{\sqrt {c+d}}\right )}{f \sqrt {c-d} \sqrt {c+d} (b c-a d)} \]

[In]

Int[Sec[e + f*x]/((a + b*Sec[e + f*x])*(c + d*Sec[e + f*x])),x]

[Out]

(2*b*ArcTanh[(Sqrt[a - b]*Tan[(e + f*x)/2])/Sqrt[a + b]])/(Sqrt[a - b]*Sqrt[a + b]*(b*c - a*d)*f) - (2*d*ArcTa
nh[(Sqrt[c - d]*Tan[(e + f*x)/2])/Sqrt[c + d]])/(Sqrt[c - d]*Sqrt[c + d]*(b*c - a*d)*f)

Rule 214

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[-a/b, 2]/a)*ArcTanh[x/Rt[-a/b, 2]], x] /; FreeQ[{a, b},
x] && NegQ[a/b]

Rule 2738

Int[((a_) + (b_.)*sin[Pi/2 + (c_.) + (d_.)*(x_)])^(-1), x_Symbol] :> With[{e = FreeFactors[Tan[(c + d*x)/2], x
]}, Dist[2*(e/d), Subst[Int[1/(a + b + (a - b)*e^2*x^2), x], x, Tan[(c + d*x)/2]/e], x]] /; FreeQ[{a, b, c, d}
, x] && NeQ[a^2 - b^2, 0]

Rule 3080

Int[((A_.) + (B_.)*sin[(e_.) + (f_.)*(x_)])/(((a_.) + (b_.)*sin[(e_.) + (f_.)*(x_)])*((c_.) + (d_.)*sin[(e_.)
+ (f_.)*(x_)])), x_Symbol] :> Dist[(A*b - a*B)/(b*c - a*d), Int[1/(a + b*Sin[e + f*x]), x], x] + Dist[(B*c - A
*d)/(b*c - a*d), Int[1/(c + d*Sin[e + f*x]), x], x] /; FreeQ[{a, b, c, d, e, f, A, B}, x] && NeQ[b*c - a*d, 0]
&& NeQ[a^2 - b^2, 0] && NeQ[c^2 - d^2, 0]

Rule 4073

Int[(csc[(e_.) + (f_.)*(x_)]*(g_.))^(p_.)*(csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_))^(m_)*(csc[(e_.) + (f_.)*(x_)]
*(d_.) + (c_))^(n_), x_Symbol] :> Dist[1/g^(m + n), Int[(g*Csc[e + f*x])^(m + n + p)*(b + a*Sin[e + f*x])^m*(d
+ c*Sin[e + f*x])^n, x], x] /; FreeQ[{a, b, c, d, e, f, g, p}, x] && NeQ[b*c - a*d, 0] && IntegerQ[m] && Inte
gerQ[n]

Rubi steps \begin{align*} \text {integral}& = \int \frac {\cos (e+f x)}{(b+a \cos (e+f x)) (d+c \cos (e+f x))} \, dx \\ & = \frac {b \int \frac {1}{b+a \cos (e+f x)} \, dx}{b c-a d}-\frac {d \int \frac {1}{d+c \cos (e+f x)} \, dx}{b c-a d} \\ & = \frac {(2 b) \text {Subst}\left (\int \frac {1}{a+b+(-a+b) x^2} \, dx,x,\tan \left (\frac {1}{2} (e+f x)\right )\right )}{(b c-a d) f}-\frac {(2 d) \text {Subst}\left (\int \frac {1}{c+d+(-c+d) x^2} \, dx,x,\tan \left (\frac {1}{2} (e+f x)\right )\right )}{(b c-a d) f} \\ & = \frac {2 b \text {arctanh}\left (\frac {\sqrt {a-b} \tan \left (\frac {1}{2} (e+f x)\right )}{\sqrt {a+b}}\right )}{\sqrt {a-b} \sqrt {a+b} (b c-a d) f}-\frac {2 d \text {arctanh}\left (\frac {\sqrt {c-d} \tan \left (\frac {1}{2} (e+f x)\right )}{\sqrt {c+d}}\right )}{\sqrt {c-d} \sqrt {c+d} (b c-a d) f} \\ \end{align*}

Mathematica [A] (verified)

Time = 0.44 (sec) , antiderivative size = 119, normalized size of antiderivative = 0.98 \[ \int \frac {\sec (e+f x)}{(a+b \sec (e+f x)) (c+d \sec (e+f x))} \, dx=-\frac {2 b \text {arctanh}\left (\frac {(-a+b) \tan \left (\frac {1}{2} (e+f x)\right )}{\sqrt {a^2-b^2}}\right )}{\sqrt {a^2-b^2} (b c-a d) f}-\frac {2 d \text {arctanh}\left (\frac {(-c+d) \tan \left (\frac {1}{2} (e+f x)\right )}{\sqrt {c^2-d^2}}\right )}{(-b c+a d) \sqrt {c^2-d^2} f} \]

[In]

Integrate[Sec[e + f*x]/((a + b*Sec[e + f*x])*(c + d*Sec[e + f*x])),x]

[Out]

(-2*b*ArcTanh[((-a + b)*Tan[(e + f*x)/2])/Sqrt[a^2 - b^2]])/(Sqrt[a^2 - b^2]*(b*c - a*d)*f) - (2*d*ArcTanh[((-
c + d)*Tan[(e + f*x)/2])/Sqrt[c^2 - d^2]])/((-(b*c) + a*d)*Sqrt[c^2 - d^2]*f)

Maple [A] (verified)

Time = 1.02 (sec) , antiderivative size = 108, normalized size of antiderivative = 0.89

method	result	size
derivativedivides	\(\frac {\frac {2 d \,\operatorname {arctanh}\left (\frac {\left (c -d \right ) \tan \left (\frac {f x}{2}+\frac {e}{2}\right )}{\sqrt {\left (c +d \right ) \left (c -d \right )}}\right )}{\left (a d -b c \right ) \sqrt {\left (c +d \right ) \left (c -d \right )}}-\frac {2 b \,\operatorname {arctanh}\left (\frac {\left (a -b \right ) \tan \left (\frac {f x}{2}+\frac {e}{2}\right )}{\sqrt {\left (a -b \right ) \left (a +b \right )}}\right )}{\left (a d -b c \right ) \sqrt {\left (a -b \right ) \left (a +b \right )}}}{f}\)	\(108\)
default	\(\frac {\frac {2 d \,\operatorname {arctanh}\left (\frac {\left (c -d \right ) \tan \left (\frac {f x}{2}+\frac {e}{2}\right )}{\sqrt {\left (c +d \right ) \left (c -d \right )}}\right )}{\left (a d -b c \right ) \sqrt {\left (c +d \right ) \left (c -d \right )}}-\frac {2 b \,\operatorname {arctanh}\left (\frac {\left (a -b \right ) \tan \left (\frac {f x}{2}+\frac {e}{2}\right )}{\sqrt {\left (a -b \right ) \left (a +b \right )}}\right )}{\left (a d -b c \right ) \sqrt {\left (a -b \right ) \left (a +b \right )}}}{f}\)	\(108\)
risch	\(\frac {b \ln \left ({\mathrm e}^{i \left (f x +e \right )}-\frac {i a^{2}-i b^{2}-b \sqrt {a^{2}-b^{2}}}{\sqrt {a^{2}-b^{2}}\, a}\right )}{\sqrt {a^{2}-b^{2}}\, \left (a d -b c \right ) f}-\frac {b \ln \left ({\mathrm e}^{i \left (f x +e \right )}+\frac {i a^{2}-i b^{2}+b \sqrt {a^{2}-b^{2}}}{\sqrt {a^{2}-b^{2}}\, a}\right )}{\sqrt {a^{2}-b^{2}}\, \left (a d -b c \right ) f}+\frac {d \ln \left ({\mathrm e}^{i \left (f x +e \right )}+\frac {i c^{2}-i d^{2}+\sqrt {c^{2}-d^{2}}\, d}{\sqrt {c^{2}-d^{2}}\, c}\right )}{\sqrt {c^{2}-d^{2}}\, \left (a d -b c \right ) f}-\frac {d \ln \left ({\mathrm e}^{i \left (f x +e \right )}+\frac {-i c^{2}+i d^{2}+\sqrt {c^{2}-d^{2}}\, d}{c \sqrt {c^{2}-d^{2}}}\right )}{\sqrt {c^{2}-d^{2}}\, \left (a d -b c \right ) f}\)	\(318\)

[In]

int(sec(f*x+e)/(a+b*sec(f*x+e))/(c+d*sec(f*x+e)),x,method=_RETURNVERBOSE)

[Out]

1/f*(2*d/(a*d-b*c)/((c+d)*(c-d))^(1/2)*arctanh((c-d)*tan(1/2*f*x+1/2*e)/((c+d)*(c-d))^(1/2))-2*b/(a*d-b*c)/((a
-b)*(a+b))^(1/2)*arctanh((a-b)*tan(1/2*f*x+1/2*e)/((a-b)*(a+b))^(1/2)))

Fricas [A] (verification not implemented)

none

Time = 1.86 (sec) , antiderivative size = 1040, normalized size of antiderivative = 8.60 \[ \int \frac {\sec (e+f x)}{(a+b \sec (e+f x)) (c+d \sec (e+f x))} \, dx=\text {Too large to display} \]

[In]

integrate(sec(f*x+e)/(a+b*sec(f*x+e))/(c+d*sec(f*x+e)),x, algorithm="fricas")

[Out]

[-1/2*((a^2 - b^2)*sqrt(c^2 - d^2)*d*log((2*c*d*cos(f*x + e) - (c^2 - 2*d^2)*cos(f*x + e)^2 + 2*sqrt(c^2 - d^2
)*(d*cos(f*x + e) + c)*sin(f*x + e) + 2*c^2 - d^2)/(c^2*cos(f*x + e)^2 + 2*c*d*cos(f*x + e) + d^2)) + (b*c^2 -
b*d^2)*sqrt(a^2 - b^2)*log((2*a*b*cos(f*x + e) - (a^2 - 2*b^2)*cos(f*x + e)^2 - 2*sqrt(a^2 - b^2)*(b*cos(f*x
+ e) + a)*sin(f*x + e) + 2*a^2 - b^2)/(a^2*cos(f*x + e)^2 + 2*a*b*cos(f*x + e) + b^2)))/(((a^2*b - b^3)*c^3 -
(a^3 - a*b^2)*c^2*d - (a^2*b - b^3)*c*d^2 + (a^3 - a*b^2)*d^3)*f), -1/2*((a^2 - b^2)*sqrt(c^2 - d^2)*d*log((2*
c*d*cos(f*x + e) - (c^2 - 2*d^2)*cos(f*x + e)^2 + 2*sqrt(c^2 - d^2)*(d*cos(f*x + e) + c)*sin(f*x + e) + 2*c^2
- d^2)/(c^2*cos(f*x + e)^2 + 2*c*d*cos(f*x + e) + d^2)) - 2*(b*c^2 - b*d^2)*sqrt(-a^2 + b^2)*arctan(-sqrt(-a^2
+ b^2)*(b*cos(f*x + e) + a)/((a^2 - b^2)*sin(f*x + e))))/(((a^2*b - b^3)*c^3 - (a^3 - a*b^2)*c^2*d - (a^2*b -
b^3)*c*d^2 + (a^3 - a*b^2)*d^3)*f), -1/2*(2*(a^2 - b^2)*sqrt(-c^2 + d^2)*d*arctan(-sqrt(-c^2 + d^2)*(d*cos(f*
x + e) + c)/((c^2 - d^2)*sin(f*x + e))) + (b*c^2 - b*d^2)*sqrt(a^2 - b^2)*log((2*a*b*cos(f*x + e) - (a^2 - 2*b
^2)*cos(f*x + e)^2 - 2*sqrt(a^2 - b^2)*(b*cos(f*x + e) + a)*sin(f*x + e) + 2*a^2 - b^2)/(a^2*cos(f*x + e)^2 +
2*a*b*cos(f*x + e) + b^2)))/(((a^2*b - b^3)*c^3 - (a^3 - a*b^2)*c^2*d - (a^2*b - b^3)*c*d^2 + (a^3 - a*b^2)*d^
3)*f), -((a^2 - b^2)*sqrt(-c^2 + d^2)*d*arctan(-sqrt(-c^2 + d^2)*(d*cos(f*x + e) + c)/((c^2 - d^2)*sin(f*x + e
))) - (b*c^2 - b*d^2)*sqrt(-a^2 + b^2)*arctan(-sqrt(-a^2 + b^2)*(b*cos(f*x + e) + a)/((a^2 - b^2)*sin(f*x + e)
)))/(((a^2*b - b^3)*c^3 - (a^3 - a*b^2)*c^2*d - (a^2*b - b^3)*c*d^2 + (a^3 - a*b^2)*d^3)*f)]

Sympy [F]

\[ \int \frac {\sec (e+f x)}{(a+b \sec (e+f x)) (c+d \sec (e+f x))} \, dx=\int \frac {\sec {\left (e + f x \right )}}{\left (a + b \sec {\left (e + f x \right )}\right ) \left (c + d \sec {\left (e + f x \right )}\right )}\, dx \]

[In]

integrate(sec(f*x+e)/(a+b*sec(f*x+e))/(c+d*sec(f*x+e)),x)

[Out]

Integral(sec(e + f*x)/((a + b*sec(e + f*x))*(c + d*sec(e + f*x))), x)

Maxima [F(-2)]

Exception generated. \[ \int \frac {\sec (e+f x)}{(a+b \sec (e+f x)) (c+d \sec (e+f x))} \, dx=\text {Exception raised: ValueError} \]

[In]

integrate(sec(f*x+e)/(a+b*sec(f*x+e))/(c+d*sec(f*x+e)),x, algorithm="maxima")

[Out]

Exception raised: ValueError >> Computation failed since Maxima requested additional constraints; using the 'a
ssume' command before evaluation *may* help (example of legal syntax is 'assume(4*c^2-4*d^2>0)', see `assume?`
for more de

Giac [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 522 vs. \(2 (103) = 206\).

Time = 0.40 (sec) , antiderivative size = 522, normalized size of antiderivative = 4.31 \[ \int \frac {\sec (e+f x)}{(a+b \sec (e+f x)) (c+d \sec (e+f x))} \, dx=\frac {\frac {{\left (\sqrt {-c^{2} + d^{2}} b {\left (c - 2 \, d\right )} {\left | c - d \right |} + \sqrt {-c^{2} + d^{2}} a d {\left | c - d \right |} + \sqrt {-c^{2} + d^{2}} {\left | -b c + a d \right |} {\left | c - d \right |}\right )} {\left (\pi \left \lfloor \frac {f x + e}{2 \, \pi } + \frac {1}{2} \right \rfloor + \arctan \left (\frac {2 \, \sqrt {\frac {1}{2}} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )}{\sqrt {-\frac {2 \, a c - 2 \, b d + \sqrt {-4 \, {\left (a c + b c + a d + b d\right )} {\left (a c - b c - a d + b d\right )} + 4 \, {\left (a c - b d\right )}^{2}}}{a c - b c - a d + b d}}}\right )\right )}}{{\left (b c - a d\right )}^{2} {\left (c^{2} - 2 \, c d + d^{2}\right )} + {\left (c^{3} - 2 \, c^{2} d + c d^{2}\right )} a {\left | -b c + a d \right |} - {\left (c^{2} d - 2 \, c d^{2} + d^{3}\right )} b {\left | -b c + a d \right |}} + \frac {{\left (\sqrt {-a^{2} + b^{2}} b c {\left | a - b \right |} + \sqrt {-a^{2} + b^{2}} {\left (a - 2 \, b\right )} d {\left | a - b \right |} - \sqrt {-a^{2} + b^{2}} {\left | -b c + a d \right |} {\left | a - b \right |}\right )} {\left (\pi \left \lfloor \frac {f x + e}{2 \, \pi } + \frac {1}{2} \right \rfloor + \arctan \left (\frac {2 \, \sqrt {\frac {1}{2}} \tan \left (\frac {1}{2} \, f x + \frac {1}{2} \, e\right )}{\sqrt {-\frac {2 \, a c - 2 \, b d - \sqrt {-4 \, {\left (a c + b c + a d + b d\right )} {\left (a c - b c - a d + b d\right )} + 4 \, {\left (a c - b d\right )}^{2}}}{a c - b c - a d + b d}}}\right )\right )}}{{\left (a^{2} - 2 \, a b + b^{2}\right )} {\left (b c - a d\right )}^{2} - {\left (a^{3} - 2 \, a^{2} b + a b^{2}\right )} c {\left | -b c + a d \right |} + {\left (a^{2} b - 2 \, a b^{2} + b^{3}\right )} d {\left | -b c + a d \right |}}}{f} \]

[In]

integrate(sec(f*x+e)/(a+b*sec(f*x+e))/(c+d*sec(f*x+e)),x, algorithm="giac")

[Out]

((sqrt(-c^2 + d^2)*b*(c - 2*d)*abs(c - d) + sqrt(-c^2 + d^2)*a*d*abs(c - d) + sqrt(-c^2 + d^2)*abs(-b*c + a*d)
*abs(c - d))*(pi*floor(1/2*(f*x + e)/pi + 1/2) + arctan(2*sqrt(1/2)*tan(1/2*f*x + 1/2*e)/sqrt(-(2*a*c - 2*b*d
+ sqrt(-4*(a*c + b*c + a*d + b*d)*(a*c - b*c - a*d + b*d) + 4*(a*c - b*d)^2))/(a*c - b*c - a*d + b*d))))/((b*c
- a*d)^2*(c^2 - 2*c*d + d^2) + (c^3 - 2*c^2*d + c*d^2)*a*abs(-b*c + a*d) - (c^2*d - 2*c*d^2 + d^3)*b*abs(-b*c
+ a*d)) + (sqrt(-a^2 + b^2)*b*c*abs(a - b) + sqrt(-a^2 + b^2)*(a - 2*b)*d*abs(a - b) - sqrt(-a^2 + b^2)*abs(-
b*c + a*d)*abs(a - b))*(pi*floor(1/2*(f*x + e)/pi + 1/2) + arctan(2*sqrt(1/2)*tan(1/2*f*x + 1/2*e)/sqrt(-(2*a*
c - 2*b*d - sqrt(-4*(a*c + b*c + a*d + b*d)*(a*c - b*c - a*d + b*d) + 4*(a*c - b*d)^2))/(a*c - b*c - a*d + b*d
))))/((a^2 - 2*a*b + b^2)*(b*c - a*d)^2 - (a^3 - 2*a^2*b + a*b^2)*c*abs(-b*c + a*d) + (a^2*b - 2*a*b^2 + b^3)*
d*abs(-b*c + a*d)))/f

Mupad [B] (verification not implemented)

Time = 15.93 (sec) , antiderivative size = 2665, normalized size of antiderivative = 22.02 \[ \int \frac {\sec (e+f x)}{(a+b \sec (e+f x)) (c+d \sec (e+f x))} \, dx=\text {Too large to display} \]

[In]

int(1/(cos(e + f*x)*(a + b/cos(e + f*x))*(c + d/cos(e + f*x))),x)

[Out]

(b*c^2*atan((b^5*c^2*tan(e/2 + (f*x)/2)*(a^2 - b^2)^(1/2)*1i - a^5*d^2*tan(e/2 + (f*x)/2)*(a^2 - b^2)^(1/2)*1i
+ b^3*d^2*tan(e/2 + (f*x)/2)*(a^2 - b^2)^(3/2)*2i + b^5*d^2*tan(e/2 + (f*x)/2)*(a^2 - b^2)^(1/2)*2i - a^2*b^3
*c^2*tan(e/2 + (f*x)/2)*(a^2 - b^2)^(1/2)*1i - a^3*b^2*c^2*tan(e/2 + (f*x)/2)*(a^2 - b^2)^(1/2)*1i - a^2*b^3*d
^2*tan(e/2 + (f*x)/2)*(a^2 - b^2)^(1/2)*3i + a^3*b^2*d^2*tan(e/2 + (f*x)/2)*(a^2 - b^2)^(1/2)*1i - b^3*c*d*tan
(e/2 + (f*x)/2)*(a^2 - b^2)^(3/2)*2i - b^5*c*d*tan(e/2 + (f*x)/2)*(a^2 - b^2)^(1/2)*2i + a*b^2*c^2*tan(e/2 + (
f*x)/2)*(a^2 - b^2)^(3/2)*2i + a*b^4*c^2*tan(e/2 + (f*x)/2)*(a^2 - b^2)^(1/2)*1i + a^4*b*d^2*tan(e/2 + (f*x)/2
)*(a^2 - b^2)^(1/2)*1i + a^2*b^3*c*d*tan(e/2 + (f*x)/2)*(a^2 - b^2)^(1/2)*2i + a^3*b^2*c*d*tan(e/2 + (f*x)/2)*
(a^2 - b^2)^(1/2)*2i - a*b^2*c*d*tan(e/2 + (f*x)/2)*(a^2 - b^2)^(3/2)*2i - a*b^4*c*d*tan(e/2 + (f*x)/2)*(a^2 -
b^2)^(1/2)*2i)/(a^6*d^2 - b^6*c^2 + 2*a^2*b^4*c^2 - a^4*b^2*c^2 + a^2*b^4*d^2 - 2*a^4*b^2*d^2))*(a^2 - b^2)^(
1/2)*2i)/(f*(a^3*d^3 - b^3*c^3 + a^2*b*c^3 - a*b^2*d^3 - a^3*c^2*d + b^3*c*d^2 + a*b^2*c^2*d - a^2*b*c*d^2)) -
(b*d^2*atan((b^5*c^2*tan(e/2 + (f*x)/2)*(a^2 - b^2)^(1/2)*1i - a^5*d^2*tan(e/2 + (f*x)/2)*(a^2 - b^2)^(1/2)*1
i + b^3*d^2*tan(e/2 + (f*x)/2)*(a^2 - b^2)^(3/2)*2i + b^5*d^2*tan(e/2 + (f*x)/2)*(a^2 - b^2)^(1/2)*2i - a^2*b^
3*c^2*tan(e/2 + (f*x)/2)*(a^2 - b^2)^(1/2)*1i - a^3*b^2*c^2*tan(e/2 + (f*x)/2)*(a^2 - b^2)^(1/2)*1i - a^2*b^3*
d^2*tan(e/2 + (f*x)/2)*(a^2 - b^2)^(1/2)*3i + a^3*b^2*d^2*tan(e/2 + (f*x)/2)*(a^2 - b^2)^(1/2)*1i - b^3*c*d*ta
n(e/2 + (f*x)/2)*(a^2 - b^2)^(3/2)*2i - b^5*c*d*tan(e/2 + (f*x)/2)*(a^2 - b^2)^(1/2)*2i + a*b^2*c^2*tan(e/2 +
(f*x)/2)*(a^2 - b^2)^(3/2)*2i + a*b^4*c^2*tan(e/2 + (f*x)/2)*(a^2 - b^2)^(1/2)*1i + a^4*b*d^2*tan(e/2 + (f*x)/
2)*(a^2 - b^2)^(1/2)*1i + a^2*b^3*c*d*tan(e/2 + (f*x)/2)*(a^2 - b^2)^(1/2)*2i + a^3*b^2*c*d*tan(e/2 + (f*x)/2)
*(a^2 - b^2)^(1/2)*2i - a*b^2*c*d*tan(e/2 + (f*x)/2)*(a^2 - b^2)^(3/2)*2i - a*b^4*c*d*tan(e/2 + (f*x)/2)*(a^2
- b^2)^(1/2)*2i)/(a^6*d^2 - b^6*c^2 + 2*a^2*b^4*c^2 - a^4*b^2*c^2 + a^2*b^4*d^2 - 2*a^4*b^2*d^2))*(a^2 - b^2)^
(1/2)*2i)/(f*(a^3*d^3 - b^3*c^3 + a^2*b*c^3 - a*b^2*d^3 - a^3*c^2*d + b^3*c*d^2 + a*b^2*c^2*d - a^2*b*c*d^2))
+ (a^2*d*atan((a^2*d^5*tan(e/2 + (f*x)/2)*(c^2 - d^2)^(1/2)*1i - b^2*c^5*tan(e/2 + (f*x)/2)*(c^2 - d^2)^(1/2)*
1i + b^2*d^3*tan(e/2 + (f*x)/2)*(c^2 - d^2)^(3/2)*2i + b^2*d^5*tan(e/2 + (f*x)/2)*(c^2 - d^2)^(1/2)*2i - a^2*c
^2*d^3*tan(e/2 + (f*x)/2)*(c^2 - d^2)^(1/2)*1i - a^2*c^3*d^2*tan(e/2 + (f*x)/2)*(c^2 - d^2)^(1/2)*1i - b^2*c^2
*d^3*tan(e/2 + (f*x)/2)*(c^2 - d^2)^(1/2)*3i + b^2*c^3*d^2*tan(e/2 + (f*x)/2)*(c^2 - d^2)^(1/2)*1i - a*b*d^3*t
an(e/2 + (f*x)/2)*(c^2 - d^2)^(3/2)*2i - a*b*d^5*tan(e/2 + (f*x)/2)*(c^2 - d^2)^(1/2)*2i + a^2*c*d^2*tan(e/2 +
(f*x)/2)*(c^2 - d^2)^(3/2)*2i + a^2*c*d^4*tan(e/2 + (f*x)/2)*(c^2 - d^2)^(1/2)*1i + b^2*c^4*d*tan(e/2 + (f*x)
/2)*(c^2 - d^2)^(1/2)*1i + a*b*c^2*d^3*tan(e/2 + (f*x)/2)*(c^2 - d^2)^(1/2)*2i + a*b*c^3*d^2*tan(e/2 + (f*x)/2
)*(c^2 - d^2)^(1/2)*2i - a*b*c*d^2*tan(e/2 + (f*x)/2)*(c^2 - d^2)^(3/2)*2i - a*b*c*d^4*tan(e/2 + (f*x)/2)*(c^2
- d^2)^(1/2)*2i)/(a^2*d^6 - b^2*c^6 - 2*a^2*c^2*d^4 + a^2*c^4*d^2 - b^2*c^2*d^4 + 2*b^2*c^4*d^2))*(c^2 - d^2)
^(1/2)*2i)/(f*(a^3*d^3 - b^3*c^3 + a^2*b*c^3 - a*b^2*d^3 - a^3*c^2*d + b^3*c*d^2 + a*b^2*c^2*d - a^2*b*c*d^2))
- (b^2*d*atan((a^2*d^5*tan(e/2 + (f*x)/2)*(c^2 - d^2)^(1/2)*1i - b^2*c^5*tan(e/2 + (f*x)/2)*(c^2 - d^2)^(1/2)
*1i + b^2*d^3*tan(e/2 + (f*x)/2)*(c^2 - d^2)^(3/2)*2i + b^2*d^5*tan(e/2 + (f*x)/2)*(c^2 - d^2)^(1/2)*2i - a^2*
c^2*d^3*tan(e/2 + (f*x)/2)*(c^2 - d^2)^(1/2)*1i - a^2*c^3*d^2*tan(e/2 + (f*x)/2)*(c^2 - d^2)^(1/2)*1i - b^2*c^
2*d^3*tan(e/2 + (f*x)/2)*(c^2 - d^2)^(1/2)*3i + b^2*c^3*d^2*tan(e/2 + (f*x)/2)*(c^2 - d^2)^(1/2)*1i - a*b*d^3*
tan(e/2 + (f*x)/2)*(c^2 - d^2)^(3/2)*2i - a*b*d^5*tan(e/2 + (f*x)/2)*(c^2 - d^2)^(1/2)*2i + a^2*c*d^2*tan(e/2
+ (f*x)/2)*(c^2 - d^2)^(3/2)*2i + a^2*c*d^4*tan(e/2 + (f*x)/2)*(c^2 - d^2)^(1/2)*1i + b^2*c^4*d*tan(e/2 + (f*x
)/2)*(c^2 - d^2)^(1/2)*1i + a*b*c^2*d^3*tan(e/2 + (f*x)/2)*(c^2 - d^2)^(1/2)*2i + a*b*c^3*d^2*tan(e/2 + (f*x)/
2)*(c^2 - d^2)^(1/2)*2i - a*b*c*d^2*tan(e/2 + (f*x)/2)*(c^2 - d^2)^(3/2)*2i - a*b*c*d^4*tan(e/2 + (f*x)/2)*(c^
2 - d^2)^(1/2)*2i)/(a^2*d^6 - b^2*c^6 - 2*a^2*c^2*d^4 + a^2*c^4*d^2 - b^2*c^2*d^4 + 2*b^2*c^4*d^2))*(c^2 - d^2
)^(1/2)*2i)/(f*(a^3*d^3 - b^3*c^3 + a^2*b*c^3 - a*b^2*d^3 - a^3*c^2*d + b^3*c*d^2 + a*b^2*c^2*d - a^2*b*c*d^2)
)